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When canonical Hamiltonians of local quantum field theories are trans- 
formed using a renormalization group procedure for effective particles, the 
resulting interaction terms are non-local. The range of their non-locality 
depends on the arbitrary parameter of scale, which characterizes the size 
of effective particles in terms of the allowed range of virtual energy changes 
caused by interactions. This article describes a generic example of the non- 
locality that characterizes light-front interaction Hamiltonian densities of 
first-order in an effective coupling constant. The same non-locality is also 
related to a relative motion wave function for a bound state of two particles. 
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1. Introduction 

Canonical quantum field theory (QFT) can be formulated using differ- 
ent forms of Hamiltonian dynamics [1]. One of the forms is distinguished 
from others by having 7 kinematical symmetries instead of usual 6. The 
distinguished form is called a front form because the Hamiltonian density 
in it is defined on a space-time hyper-plane that is swept by the wave front 
of a plane wave of light. The hyper-plane is called light front (LF). The 
7th symmetry is invariance with respect to the Lorentz boosts along the 
direction of motion of the wave. The additional symmetry has many conse- 
quences. For example, one can describe the relative motion of constituents 
of a bound state for arbitrary motion of the bound state as a whole, achiev- 
ing a connection between the rest frame image of the system with its image 
in the infinite momentum frame, and the ground state (vacuum) problem 
in a theory is posed in a new way because there is no spontaneous creation 
of particles from empty space [2] . A canonical formulation of the standard 
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model (SM) using LF hyperplane for quantization of fields can be found 
in [3]. 

Canonical Hamiltonians of local theories are singular operators. They 
require regularization and renormalization. In particular, a renormalization 
group (RG) procedure produces Hamiltonians that are called renormalized. 
Renormalized Hamiltonians depend on a RG scale parameter, here called 
A, and are no longer local. The size of their non-locality corresponds to the 
size of A. A renormalized Hamiltonian of some scale A can also be called 
effective up to the scale A, since it is equivalent to the Hamiltonian of the 
original theory and provides an optimal setup for calculating observables 
that concern phenomena whose description does not require a resolution 
of effects beyond the resolution implied by the size of A. Effective Hamil- 
tonians do not depend on regularization but they do depend on the scale 
A. This article discusses lowest-order non- localities that characterize effec- 
tive Hamiltonian densities on the LF. Besides general interest in effective 
non-local interaction densities in quantum field theory, LF non-localities are 
intriguing because of their necessarily relativistic nature. 

Consider the case of QCD, which is a part of the SM. The LF power- 
counting for Hamiltonian densities in QCD indicates [4] that a large number 
of operators may contribute to a renormalized Hamiltonian. The power- 
counting leads to a complex set of operators because one works in the 
Minkowski space, the direction of z-axis (the direction of motion of the plane 
wave of light that defines the LF is conventionally chosen to be against z- 
axis) and the directions transverse to z-axis are treated differently, and, ab 
initio, one has to deal with operators that have important matrix elements 
arbitrarily far off energy shell. In summary, effective LF Hamiltonians are 
not expected to have a structure that can be guessed easily. One needs a 
method to derive them. 

The method discussed here is the renormalization group procedure for 
effective particles [5] . The name of the method is abbreviated RGPEP. The 
method assumes that finite parts of counterterms can be fixed using predic- 
tions for observables that follow from the calculated effective Hamiltonians. 
The observables may include properties of bound states even if a procedure 
for calculating the Hamiltonians is itself carried out only in an order-by- 
order expansion in powers of an effective coupling constant. This can be 
seen on the example of a Coulomb interaction in the Schrodinger Hamilto- 
nian for an electron and a proton. The interaction is proportional to the 
first power of a and yet it predicts properties of the electron-proton bound 
states whose wave functions have no expansion in powers of a around 0. 
Thus, one may hope to learn about non-perturbative solutions of a theory 
by solving eigenvalue problems for Hamiltonians derived using RGPEP even 
in low orders of perturbation theory. 
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In the ultraviolet regime of QCD, one obtains non-local vertices that ex- 
hibit asymptotic freedom [6, 7] in the dependence of the effective coupling 
constant in them, g\, on the running RG parameter A [8]. At the same 
time, the parameter A limits from above allowed energy changes in the in- 
teraction vertices 1 . This means, through the connection between energy and 
momentum and uncertainty principle for momentum and position variables, 
that the vertices resulting from RGPEP are non-local and the range of the 
non-locality is inversely related to the size of A. In the ultraviolet regime, 
one has A 3> Aqcd, where Aqcd is the specific momentum scale parameter 
that physically characterizes QCD in the scheme of RGPEP. The limit of 
infinite A, if obtained, could define a local theory. 

At the opposite end of the A scale, there lies the infrared regime charac- 
terized by A ~ Aqcd- In the infrared regime, one has to deal with Hamilto- 
nian terms that are involved in the formation of bound states. These terms 
are new in the sense that they may significantly differ in appearance from 
terms implied by a classical Lagrangian in a canonical definition of QCD. 
In this regime, RGPEP is designed to guarantee that only invariant mass 
changes up to A ~ Aqcd are allowed to occur (see below). This guaran- 
tee is associated with the characteristic non-localities of RGPEP that are 
described in this article. 

In the infrared regime, one also expects that the masses of effective par- 
ticles receive dynamic contributions order Aqcd- Therefore, the effective 
dynamics is expected to be limited to slow relative motion of interacting ef- 
fective particles of masses presumably not smaller in size than about Aqcd- 

The infrared regime is not precisely understood. Nevertheless, on the 
basis of phenomenology of hadrons and strong interactions, one expects 
formation of constituent quarks and gluons in this regime and RGPEP is 
designed to help in identifying the mechanism that could lead to their for- 
mation as effective particles in QCD. 

RGPEP equations would have to be solved non-perturbatively in order 
to provide an exact effective Hamiltonian of QCD for any A. Exact cal- 
culations are hardly possible and it is not shown that the eigenvalues of 
LF Hamiltonians, equivalent to masses squared of physical states, must be 
positive. A negative mass squared eigenvalue would have to be explained. 
There is no explanation other than a failure of a theory or a method applied 
to seek solutions, such as a wrong identification of a theory ground state. 
However, since the LF RGPEP allows us to calculate effective Hamiltoni- 
ans order-by-order in powers of g\ without making assumptions about the 
QCD ground state [4], such as the assumption of vacuum condensates [9], 

1 In the LF form of dynamics, A of RGPEP limits changes in the invariant mass of 
interacting particles instead of changes in their energy, and only a subset of particles 
that are involved in a single action of an effective Hamiltonian counts in the difference. 
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the structure of the LF non-local interactions that result from RGPEP in 
the order-by-order calculations are of interest as elements of the unknown 
realm. This article provides a description of the leading non-localities that 
show up already in the first-order terms, independently of additional non- 
local effects that build up on top of the leading ones in higher-order terms. 

Interest in the non-locality that characterize RGPEP is also motivated 
by a general observation that the non-locality contributes to a replacement 
of an abstract local theory (the one that requires regularization) by a per- 
haps more realistic effective one, whose scale parameter A can be adjusted 
in order to obtain a Hamiltonian whose physically most important terms 
look simplest possible [10]. In this sense, the non-locality is obtained on the 
basis of RG ideas with no need for invoking any new degrees of freedom, 
such as, for example, the ones introduced in the case of a non-local string 
picture of particles 2 or other forms of substructure, perhaps even includ- 
ing additional symmetries. From this point of view, the non-local structure 
of interactions that RGPEP may lead to is of more general interest than 
QCD alone. In fact, the first-order non-locality discussed in this article oc- 
curs in QFT in the same form irrespective of many details concerning spin 
and other quantum numbers of effective particles, such as isospin, flavor, or 
color. For example, it would be the same in perturbative quantum gravity. 

Section 2 reviews qualitative features of non-local interaction vertices 
in lowest-order RGPEP in momentum space. The main Section 3 discusses 
non-locality in space-time in several subsections, starting in Section 3.1 from 
the question of how non-local interactions become local when A — > oo. The 
non- locality for finite A and particle mass m — > is discussed in Section 3.2. 
Non-locality for A ~ m is described in Section 3.3. A comparison with a 
non-relativistic theory is provided in Section 3.4. Section 3.5 explains the 
relationship between the relativistic non-localities and 2-body bound-state 
wave functions. Section 4 concludes the paper with a summary and a few 
comments. An Appendix in three parts provides some details of calculations. 

2. Non-local vertices in momentum space 

Let a regularized Hamiltonian with counterterms for some QFT be de- 
noted by H and let a denote creation or annihilation operators of bare 
particles in this Hamiltonian. For example, the bare particles in the case of 
LF Hamiltonian of QCD are quanta of canonical quark and gluon fields in 

2 Questions concerning gravity and its relation to quantum mechanics require a sep- 
arate discussion. Also, conditions of causality may constrain Hamiltonian non- 
localities and it is not known if RGPEP can satisfy them automatically in gauge 
theories. The author would like to thank J. Lukierski for a comment concerning the 
latter issue. 



SDGlazekNonLocalityEPS printed on June 22, 2010 



■5 



the gauge A + = A + A 3 , with some arbitrarily large cutoffs on momenta. 
The cutoffs are meant to be arbitrarily large in the sense of tending to the 
limit of removing the regularization. 

In RGPEP, one introduces effective particles of scale A through a trans- 
formation U\ (see below) whose construction involves the entire H [5]. 
Counterterms introduced in H allow one to construct a formal expansion 
in a series of powers of the effective coupling constant g\, but it should 
be stressed that even the Hamiltonian of order is not fully known yet 
in important theories [8, 11]. Since the counterterms are constructed us- 
ing RGPEP (to identify their structure) and results for observables (to fix 
the unknown finite parts), the complete transformation U\ can be made 
unitary in an order-by-order construction if one can simultaneously change 
the initial H by including counterterms also order-by-order. However, some 
characteristic terms can be introduced partly in a non- perturb at ive way. 

The key example is provided by a quark mass-squared term to which 
one can add a term proportional to Aqcd ~ A 2 exp (— b/a\), where b is a 
constant and a\ = g\/(Ait). Such addition contributes to an expansion in 
powers of g\ around g\ = but the calculation of U\ can be carried out in 
RGPEP for arbitrary values of particle masses. Then, on the one hand, mass 
terms can be adjusted by matching theoretical predictions with observables. 
On the other hand, the question of how a non-perturbative dynamics de- 
scribed by the calculated Hamiltonian relates mass terms to observables re- 
quires investigation. Other examples of inclusion of non-perturbative effects 
in a LF Hamiltonian density are provided by the arrangement of couplings 
according to the rules of coupling coherence that keeps track of symmetries 
[12], including the arrangement of operators allowed by power-counting so 
that the resulting set of terms corresponds to a theory with a spontaneously 
broken symmetry [4]. 

The example of a mass term in Hq is important in discussion of non- 
local Hamiltonian densities because eigenvalues of Hq are used to define 
the RGPEP factors responsible for the non-locality. The non-locality of 
effective QCD depends on the quark and gluon mass parameters. Since 
QCD promises to generate contributions to the effective particle masses 
from Aqcd, it should be stressed that the parameter Aqcd primarily char- 
acterizes perturbative A-dependence of effective theories for large A. In 
theories with large A, small mass terms may be treated as negligible. In 
this context, the variation of non-locality of effective LF QCD Hamiltoni- 
ans over a large range of A appears related to the question of generation 
of particle masses in local theories with formal symmetries considered valid 
for strictly massless particles (chiral symmetry). This article discusses non- 
localities for different ratios of a mass parameter m to the scale parameter 
A, irrespective of the value of parameters like Aqcd- 
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The RGPEP operation U x mentioned above, transforms creation and 
annihilation operators by a rotation, 

ax = U x aU{. (1) 

The corresponding Hamiltonian operator H x is constructed to be the same 
as H, 

H x (a x ) = H(a). (2) 

Consequently, H\(a\) obtainable in perturbative RGPEP is assumed to 
be a combination of products of operators a\ with coefficients c\ that are 
different from coefficients c of corresponding products of operators a in 
H(a). RGPEP provides differential (or algebraic) equations that produce 
expressions for the coefficients c x in H x [5]. Namely, from 



and the condition Uoo = 



H x (a) = U{H(a)U x , 
1, one obtains 

-^H x (a) = [T x ,H x (a)}, 
#oo(a) = H(a), 



(3) 

(4) 
(5) 



where 

n = -vl£n- («) 

Therefore, the evolution of coefficients c x with A is determined by T x ; see 
Eqs. (2.28) and (2.29) in [5]. 

Since T x vanishes when interactions vanish, the Hamiltonian H x (a) can 
be expanded in powers of the interaction strength. Suppose the interaction 
strength is parameterized by a suitably chosen coupling constant g x . In the 
case of QCD, the limit of vanishing g x for any fixed value of A can also 
be seen as a limit of vanishing Aqcd- In this limit, the non-perturbative 
terms mentioned above, proportional to positive powers of Aqcd, are also 
vanishing. 

The first-order terms in H x (a x ), i.e., terms proportional to the first 
power of g x , are functions of operators a x . As such, they have the same 
form as the canonical interaction terms proportional to the bare coupling 
constant g have as functions of the b are operators tt, containing Cqo — c. 
The only difference is that g is replaced by g x and c is replaced by c x = f x c, 
where f x is a vertex form factor of RGPEP [5]. Precisely this form factor 
introduces the nondocality studied in this article. 
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The first step in studying the non-locality resulting from the form factor 
fx is to define a generic example of the interaction term in which fx appears. 
For this purpose, one can observe that in physically important canonical 
local theories, such as gauge theories or Yukawa theory, the Lagrangian 
interaction densities of first-order in the coupling constants contain products 
of three fields evaluated at the same space-time point. For example, fermion 
fields tp are coupled with gauge boson fields A in a product of the type 
g'tpfii/j, and they are coupled with scalar fields 4> through the product of the 
type gtjxjnf). Non-Abelian gauge fields are coupled to themselves through a 
product of the type gTr d^AJ^A^ , A v \. Therefore, for the present discussion 
of first-order non-localities, it is sufficient to consider space-time operator 
densities that are products of three fields. 

Another observation is that fx only depends on the change of energy 
across an interaction Hamiltonian Hxi = H\ — Hxo- More precisely, in LF 
dynamics, f\ depends on the change of an invariant mass of the interacting 
particles across the interaction. 3 This means that the first-order non-locality 
structure does not depend on spin, isospin, flavor or color variables. In other 
words, the first-order RGPEP evolution of coefficients cx with A is limited 
to variation of the range of allowed changes in the invariant mass and this 
change depends only on the momenta and masses of the interacting particles. 

According to these two observations, generic features of the first-order 
non- locality can be studied in the case of a product of three scalar fields. 
Conclusions regarding first-order non-locality in interactions of more com- 
plex fields will be the same as for scalars, except for additional algebraic 
factors or derivatives that originate directly from the vertices of correspond- 
ing canonical theories. 

Consider the classical canonical LF interaction Hamiltonian for a real 
(chargeless) scalar field tp( x ): 



where x~ = x° — x 3 , x 1 - = (x 1 ,^ 2 ), and the LF is defined by the condition 
x + = x° + x 3 = 0. Hermitian quantum field ip(x) is composed of creation 
and annihilation operators for bare particles, 



3 At the same time, three components of a total momentum, P + — P° + P 3 and 
P x = (P 1 , P 2 ), are preserved and the invariant-mass change is invariant with respect 
to 7 Poincare transformations that preserve the LF hyperplane. 
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where 

a- P = al p , (10) 

cfip = dp + d 2 p ± , and all three integrals over momentum variables extend 
from — oo to +oo. Single commutation relation 

[a p ,a q ] = 2p+{2it) 3 5 3 {p + q), (11) 

contains three corresponding commutation relations: one for two creation 
operators when both p + and q + are negative, one for two annihilation oper- 
ators when both p + and q + are positive, and one for one creation operator 
and one annihilation operator, when p + and q + have opposite signs. The 
commutation relation implies 

[iP(x),d + iP(y)} = iS 3 (x-y), (12) 

where d + = 2d/dx~. One also has 

a p = \p + \ j d 3 xe +ipx i;(x), (13) 

where d 3 x = dx~d 2 x ± and all three integrals over position variables extend 
from — oo to +oo on the LF. 

The above notation differs from the standard one [13, 14]. The difference 
is that the creation and annihilation operators are distinguished solely by 
the sign of p + in a p and the integration over momentum component p + is 
not limited to only positive values. The colon sign in Eq. (7) denotes normal 
ordering. The normal ordering is defined using Feynman's convention [15] 
with the ordering parameter set equal to p + that ranges from — oo to +oo. 
In this convention, it is understood that an operator a Pl stands to the left 
of the operator a P2 when pf < p^. Otherwise, the order is reversed. Thus, 
all annihilation operators a p , which by definition have p + > 0, stand to the 
right of all creation operators, which by definition are a p with p + < 0. 

In terms of operators a, the local, canonical interaction Hamiltonian 
reads 

Hi = g J[piP2P3] 2(27r) 3 5 3 (p 1 + p 2 +P3) ■ 

a Pl a p2 a p:i : . (14) 

From this expression, RGPEP produces an effective interaction Hamiltonian 
of first order in the form [5] 

H\i = g\ j [P1P2P3] 2(27r) 3 5 3 (pi + p 2 + P3) fx ■ (15) 
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where 

fx = e -(AA4VA 2 ) 2 (16) 

is the source of non-locality of the effective vertex. The argument of the 
form factor f\ is the difference between the invariant masses of annihilated 
particles, M 2 , and created particles, M 2 . Namely, 

AM 2 =M 2 a -M 2 c , M 2 a = P 2 , M 2 C =P 2 , (17) 

where the total momentum four- vectors for annihilated and created particles 
are 

Pa = fl &)Pii P c = Y. (~Pt)Pi, Pi = V ^\ m \ (18) 
i=l i=l Pi 

If a Hamiltonian density term contained n > 3 fields, the sums over mo- 
menta in the above expressions that define fx as function of a change in the 
invariant mass of interacting particles would extend up to n instead of only 
3, with no other change. 

3. Non-local vertices in position space 

In the first step of defining effective LF Hamiltonian densities on the 
LF, we introduce effective quantum field operators. This is done in analogy 
with Eqs. (8) to (13) for bare fields. Namely, 

= J[p]ax P e- tpx = 1>{(x), (19) 

where 

ax- P = a[ p , (20) 

[a Xp ,ax q ] = 2p + (2irf5 3 (p + q), (21) 

ax P = \p + \ J d 3 xe^^ x (x). (22) 

Thus, the effective field commutation relations remain the same irrespective 
of the value of A, 

lMx),d + My)} = i^{x-y). (23) 
Using Eq. (22), the effective interaction Hxi is obtained in the form 

Hxi = 9xJ [P1P2P3] 2(2tt) 3 5 3 ( Pl +p 2 + P3 ) fx \ptP2Pt\ (24) 
x J d 3 x 1 d 3 x 2 d 3 x 3 e +i( P lxl+P2X2+P3X ^ :Mxi)^\{x2)Mxz) ■ ■ 
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This means that 

H\i = g\ j d 3 x 1 d 3 x 2 d 3 x 3 f x (x 1 ,x 2 ,x 3 ) : ipx(xi) ^\{x 2 ) ^\{x 3 ) : ,(25) 

where the space-time non-locality of the effective interaction Hamiltonian 
density on the LF is described by the function 

f x (x 1 ,x 2 ,x 3 ) = y [pip 2 p 3 ] 2(2vr) 3 S 3 (p 1 +p 2 +p 3 )\pipipi\ 

Evaluation of this integral leads to the main results concerning non-locality 
of effective LF Hamiltonians. But in order to introduce relevant concepts we 
first discuss the issue of how the non-local interactions become local when 
A -> 0. 

3.1. How non-local interactions become local when A — > oo 

Locality of an interaction Hamiltonian in the limit of A — > oo is here 
understood as the following feature: matrix elements of the interaction 
Hamiltonian vanish between states of effective particles corresponding to 
RGPEP scale A if the wave functions of these particles in these states have 
supports separated by a distance r that is not smaller than some fixed but 
arbitrarily small distance tq so that rX — > oo. In other words, locality in 
the limit A — > oo means that matrix elements between all states that are 
experimentally separable in space vanish and all experimentally accessible 
separations r are considered different from by no less than some very small 
but fixed amount vq when A — > oo. Every local theory may be discovered 
inapplicable in physics when vq is reduced below certain value. In such case, 
ro designates the limit of physical applicability of that local theory. 

We begin by showing how Eq. (26) produces an interaction that becomes 
local when A — > oo. Formally, the notation introduced in Eqs. (9), (19), and 
(22) produces this result pointwise in an obvious way: limA->oo fx = /oo = 1, 
factors l/|p + | in the integration measures are cancelled by the factors of \p + \ 
in the expressions for a p = a^p = limA-s.oo a\p, and the remaining integrals 
produce 

foo(x 1 ,x 2 ,x 3 ) = S 3 (x 1 - x 3 ) 5 3 (x 2 - x 3 ) . (27) 

This result leads through Eq. (25) to Eq. (7), with operators tpoo(x) = 
limA->oo ipx(x) = ip(x) as a consequence of a = Oqo = limA->oo o,\. It is 
also understood that limA->oo 9\ = 9oo and differs from g only by the 
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value implied by a coupling constant counterterm. In this reasoning, the 
space-time variable X3 in Eq. (25) plays the role of variable x in Eq. (7). 

On the other hand, this mechanism appears to require further expla- 
nation because particle momenta p + are limited to only positive values in 
the standard LF notation for canonical quantum fields ip(x). The same 
standard notation could be used in the case of effective fields if)\ (x) , and all 
integrals over p + would only extend from to +00. Instead, in our notation, 
the momentum conservation 5-function in Eq. (15) forces one of the mo- 
menta pi, P2 and P3 to have an opposite + -component to the two others 4 . 
According to Eq. (22), this means that the only terms that contribute are 
those in which one particle of positive p + is annihilated and two particles 
of positive p + are created, or two are annihilated and one is created. There 
are 3 terms of each type. Therefore, the effective interaction equals 



f = \Jo 2pT(2ir)J 



d 2 P j 



(2vr) 2 

2(2vr) 3 S 3 ( Pl + P2 - p 3 ) f x (a[ pi a[ p2 a Xp3 + h.c) , (28) 



where 

f x = e -(>tW) 2 M 4 . (29) 



The sign of normal ordering is not needed. Thus, one can write 



H 



XI 



/ 



"ft r^-i 

l\Jo 2(2vr)i 



d 2 pj 

(27T)* 



2(27r) 3 «5 3 (p 1 +p 2 -p 3 ) fx 



d 3 xid 3 x 2 d 3 X3 



e -i(p 1Xl+P2X2 -p 3X3 ) ^ A ( Xl )^ A ( X2 )^ A ( X3 ) + h.c. 



(30) 



without normal ordering. It is visible that the momentum components pj~ 
and P2 are all always positive and they sum up to an always positive p%. 
This implies that one can write pf = zp% , p\ = (1 — z)pf , and the range of 

4 The case that all three momenta have components p + — is excluded by regulariza- 
tion in a canonical theory by demanding that \p + \ in the Fourier expansion of every 
field in a Hamiltonian density is greater than some positive infinitesimal constant 
e + . The strong limit of e + — > is immediately taken in all terms obtained from the 
integration on the LF in the theory regularized with an infinitesimally small e + . Sub- 
sequent transformation U\ replaces operators a with a\ preserving their momentum 
labels. 
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integration over z is from to 1. Moreover, introducing parameterization 

pf = zP+, pi = zP^ + q^, V ~ x = P i 2 +f , (31) 

pt = (l-z)P\ pi = (l-z)P ± - q \ p~ 2 = ll_ + ™p + , (32) 

which is a standard way of parameterizing relative motion of interacting 
particles 1 and 2 with variables z and q 1 - in canonical LF dynamics, one 
obtains the invariant mass squared of particles 1 and 2 in the argument of 
the form factor fx equal 

(33) 

The above relations hold independently of the value of A. Changing inte- 
gration variables to P + , P- 1 , z, and q- 1 -, one obtains for x + = that 

f°° dP+ f d 2 P ± f 1 dzP+ f d 2 q ± f 
Hxi ~ 39X J 2(2^)7 Wf], 2(27) J Wf h 

x J d 3 Xl d 3 x 2 d 3 x 3 [e- iX ij x (xi)iJ\(x2)ipx(x 3 ) + h.c] , (34) 

X = (zP + q)x 1 + [(l-z)P-q]x 2 -Px 3 . (35) 

But this result means that 

Hxi = 9\ j d 3 Xl d 3 x 2 d 3 x 3 fx(x 1 ,x 2 ,x 3 ) ip\(xi)ip\ {x 2 )i)x{x 3 ) + h.c. , 

(36) 



where 



j, , [°° dP+ [ d 2 P ± f 1 dzP+ f d 2 q L , 

/A(X1 ' X2 ' X3) = 3 7o 2wl Wfh 



e- jX . (37) 

There is also no need for the sign of normal ordering in Eq. (36), since the 
signs of momentum variables in the integration automatically put creation 
operators to the left of annihilation operators. 

Eqs. (36) and (37) should be compared with Eqs. (25) and (26). The 
issue to clarify is how the integration over z only from to 1 and over 
only positive values of P + leads to a local interaction when A — > oo. The 
limited range of integration in momentum space seems to always require 
some smearing of interaction in position space and the variable z is limited 
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in Eq. (37) to the range from to 1, irrespective of the size of A. How- 
does locality emerge in the limit of A — > oo? Details of relevant reasoning 
are collected in Appendix A. Here we only provide a description of the 
connection between the support of fx as function of z between and 1 and 
the locality of H\j in the limit A — > oo. 

The non-locality of interaction Hamiltonian in Eq. (25) is described by 
the function fx(xi,X2,Xs) defined in Eq. (26) by an integral over momentum 
variables. Eq. (A. 12) in Appendix A shows that 

d 3p y+oo d ^\ P +\ y d 2 K ± _^ 



A(*i,* 2 ,*3) = J ^ J _ Wj / ^ , (38) 

fx - A(C,«) =e-W/*r, (39) 



\AM 



K ±2 + m 2 2 ] 1 + |C| + |1-C| 



LC(i-C) '"J 2 • ^ 

Y = ((P + k)xi + [(1 - QP ~ k]x 2 - Px 3 . (41) 

The Hamiltonian obtained by integrating a normal-ordered product of three 
quantum fields ip(xi), ip(x2), ip(xs) with function f\(xi, X2, x%) over LF 
hyperplane in Eq. (25), is the same as the Hamiltonian obtained in Eq. 
(36). In Eq. (38) for f\(xi, X2, £3), integration over the range of negative 
P + corresponds to the sign of Hermitian conjugation, /i.e., in Eq. (36). 
Integrations over ( < and ( > 1 contribute the same operator as the 
integration over £ between and 1 does. Thus, the three regions contribute 
the factor 3 in Eq. (34). This means that the factor 3 appears in front of 
one and the same operator that is obtained in turn from each of the three 
regions equally. Appendix A describes changes of variables that exhibit the 
equivalence of these three regions of integration over (. 

Fig. 1 shows the function fx = fxiCn- 1 ) = exp [-(AA4 2 /A 2 ) 2 ], where 
AA4 2 as function of £ and is given in Eq. (40), for six different values of 
A: 50, 10, 4, 3, 2, and 1.5, in units of the mass m. Momentum k l is also given 
in units of m. The six three-dimensional plots are approximate because of 
a limited resolution of figure drawing. In particular, fx = for £ = and 
( = 1 for all values of A, except that the region where fx is very small 
decreases in size and relevance in the Fourier integral when A — > 00. This 
means that the non-locality of the effective interaction Hamiltonian H\i can 
be systematically considered only under the assumption that the domain 
of Hxi does not contain states with wave functions that are significantly 
singular as functions of £ at and 1. We assume that the domain of Hxi 
obeys this condition, on the basis of our expectation that a non-zero mass- 
term leads to a spectrum that satisfies this condition. 

The middle bumps on the six plots in Fig. 1 show the RGPEP form 
factor fx in the range < £ < 1. In this range, ( is the same as the 
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Fig. 1. The RGPEP vertex form factor fx = f\((, k) that is integrated over mo- 
menta in Eqs. (26) and (38) to produce the non-local interaction density in H\j 
in Eq. (25), for six values of A. m is the mass of interacting particles. The three 
bumps correspond to the regions C<0,0<£<1, and ( > 1, all three contributing 
equally to H\i. When A —¥ oo, f\ tends pointwise to 1. Limited resolution of figure 
drawing misses the fact that /a(0, k) = f\(l, k) = for m > (see the text). 



parameter x typically used in LF notation for + -momentum fractions. The 
other two bumps that are visible in all six plots in the regions £ < and 
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C > 1, combine together with the bump in the middle to form a function 
that approaches 1 pointwise when A — > oo. This is how one obtains a local 
interaction in this limit. 

It is visible in Fig. 1 that for A — > oo one obtains interaction whose 
strength approaches a constant in momentum space, and hence tends to a 
point-like interaction in position space. When the position probe one uses 
has a resolution considerably below 1/A, the non-locality corresponding to 
A remains invisible. The greater A, the smaller distance scales at which the 
interaction still appears to be local, i.e., the interaction density spreads only 
over regions that are smaller than the probe can resolve. 

It is also visible in Fig. 1 that the smaller A the greater overall sup- 
pression of the strength of the interaction by the form factor f\. This 
suppression can be reinforced (triviality) or reduced (asymptotic freedom) 
by the variation of the coupling constant g\ when A decreases. 

3.2. Non-locality for m — > 

High-energy dynamics involves interactions in which spatial momenta 
of particles are typically very large in comparison with their masses. In 
such circumstances, it may be useful to neglect the masses. One can also 
consider masses that are negligible in comparison with spatial momenta 
of interacting particles, even if the latter are not very large on the scale 
of momenta measurable in laboratory. This section discusses non-locality 
of first-order RGPEP three-particle vertices in the case when masses are 
neglected. 

The case of m = can be seen as resulting from the limit m — > in 
the sense that the range of relative momenta accessible to particles created 
or annihilated in a single act of interaction is much greater in size than m. 
Since the invariant mass depends on q ±2 + m 2 , large momentum q ± makes 
a small mass parameter m irrelevant for the value of the form factor f\. 
But for small momentum q^, the mass m may be important because m 2 is 
divided by z or 1 — z in the invariant mass. For q 1 - = 0, this means that an 
arbitrarily small m matters in the value of f\ wherever m 2 /z or m 2 / (1 — z) 
are not negligible in comparison with A. For m <C A, this happens only for 
extreme values of z. Therefore, one can expect that the limit m — > is 
equivalent to the result of setting m = with the exception of contributions 
from the end points in z while \q L \ < m. At the end points, the case m = 
qualitatively differs from all cases with m > 0. The former case allows for 
non-zero contributions from the end points when q 1 - — > 0, and all the other 
cases do not. 

One can gain some understanding of the non-locality corresponding 
to A > m by neglecting m entirely and replacing the form factor fx = 
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exp[-(A^i2/A) 4 ] by fx = exp[-(.Mi2/A) 2 ] in Eq. (37). Such replacement 
preserves a non-local nature of the vertex qualitatively and greatly simplifies 
calculations. The simplified case is discussed in this section. 

Under the simplifying assumptions, the non-local vertex function in Eq. 
(37) reads 

j , , fd 3 P9(P + ) f 1 dzP+ f d 2 q A 

M*> x * x *) = 3 J -ip-/ WiN 



e \*z(i-z) 

x e -iP[zx 1 + (l-z)x 2 -x 3 ]+iq ± (x^-x^) ^ rgp. 

Eq. (B.8) in Appendix B provides the result of integration over momenta 
q ± , P ± , and P + . 5 The vertex function is invariant under translations and 
depends on 6 relative position co-ordinates. For example, when one identi- 
fies a point of reference with the argument of the field labeled 2, the result 
takes the form 

f\(xi,x 2 ,x 3 ) = A 6 5-[A(xi - x 2 ), \{x 3 - x 2 )] , (43) 

where 



/' 

Jo 



f*P-T) - T^l d*z(l-*, 5 2 {T ±_ zp ± )e - Z (l- Z ) P ^ (J|) 

16ir 6 J (t - zp + ie) 1 



The function g{p,r) is invariant with respect to rotations around z-axis. It 
is different from only for being a fraction of p^. 

In order to visualize the non-local vertex function f\(x\, x 2 , x%), it is 
sufficient to consider it a function of r on the LF hyper-plane of p. More- 
over, since t 1 - = tp 1 - with t = |t- l |/|/t l | in the range < t < 1, one can 
visualize f\(x±, x 2 , X3) by drawing it on the two-dimensional plane of only 
two variables: and t~ . Namely, using the decomposition 

T L = tp ± + sn ± , (45) 

with n ± p ± = and n^n 1 - = 1, one can write 

9(P,r) = 5(8)^h„{t,T-), (46) 



where 



(r~ - tp~ + ie) 
< t = \r L \/\p ± \ < 1. (48) 



M*> O = —- i- ■% —JTTi » ( 47 ) 



Integration over P+ requires a regularization that results in the presence of ie in 
Eq. (44); see Appendix B. Despite the presence of ie, the resulting interaction 
Hamiltonian is Hermitian because Eq. (36) contains a term with f\ and a conjugated 
term, denoted by /i.e. 
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Fig. 2. Modulus of the non-local interaction density for massless particles, 
f\(xi,X2,Xs) in Eq. (42), drawn as function of Xx% and Xx% on the space-time 
plane containing x^ — x^ in the LF hyper-plane defined by the condition x + = 0, 
in terms of the function h p in Eq. (47) multiplied by \p | for visualization of the 
end-point enhancement, for Xx2 — and several choices of the points Xxi with 
Xxi = X\x± | = 1, 4, 6, and 10, as indicated. The end-point enhancement for 
massless particles when X3 approaches x\ or X2 shows up when the dimensionless 
distance X\x± — x^\ exceeds 4. For further explanation, see Eqs. (43) to (48) and 
the text. 



For every choice of p = X(x± — X2) = (p ,p~), the function h p (t,T~) is a 
function of < t < 1 and t~ . Figs. 2 and 3 show examples that illustrate 
generic features of f\(x\, X2, X3). 

Fig. 2 illustrates how the non-local interaction strength on the LF de- 
pends on the position of the point X3, where particle number 3 is annihilated 
(created), for a given distance between the two points x\ and X2, where par- 
ticles 1 and 2 are created (annihilated). The distances are measured in di- 
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mensionless units that result from multiplication of position coordinates by 
the RGPEP scale parameter A. To obtain the Hamiltonian corresponding to 
A, the density h p whose modulus times \p^\ is shown in Fig. 2, multiplied by 
the factors present in Eq. (46) to obtain f\(x±, X2, x^), is integrated with 
the product of three effective fields ip\(x) using f \ 3 d 3 x\ f X 3 d 3 X2 f d 3 x%, 
seeEq. (36). 

For the purpose of drawing, Fig. 2 is artificially modified in so far that 
the infinitesimal parameter e — > 0+ is set to 1/5. This substitution causes 
that the rise of the modulus of f\(xi, X2, X3) due to the square of x^ — x% — 
(xi — x^ ) I x% — %2 I / 1 xj~ — X2 I crossing in denominator is limited to 1 /e 2 = 
25, instead of reaching 0+ 2 . The artificial modification preserves generic 
features of f\{x\, X2,Xs)- Namely: a) f\(x±, X2, £3) pointwise decreases with 
distance between xj- and x% , b) it quickly vanishes outside the region where 
X3 is near a straight line connecting x\ with X2 on the LF, and c) it is notably 
spread toward the end-points when \\x\; — x^\ exceeds 4. The number 4 
results from the factor e - * 1 - 1- *^ 7/4 in Eq. (47). This factor is varying 
slowly until \p^\ 2 exceeds the inverse of maximal value of t(l — t)/4, which 
is 16. Thus, when \p^~\ exceeds 4, the density f\{x\, X2, X3) begins to be 
suppressed in the middle between x\ and X2- When \p ± \ increases far above 
4, the interaction density favors X3 near x\ or X2, being squeezed in the 
middle between x\ and X2- 

Fig. 2 illustrates the feature that X3 must lie near a line connecting x\ 
and X2 on the LF in the cases where |x7/ — x^\ ~ \x± — x^\- Fig. 3 shows 
that the same happens in other cases, by providing examples for different 
orientations of x\ — %2 on the LF. Again, e is set to 1/5 in order to avoid 
infinite values on the line connecting points x\ and X2, which is the same 
trick of convenience in drawing that was used earlier in Fig. 2 and is further 
explained below. 

Both Figs. 2 and 3 display partly jagged shapes of the drawn functions, 
which requires explanation. For this purpose, consider the function 

f{x) = -( x + ie )- 2 , (49) 

in the vicinity of x = 0. Fig. 4 displays the modulus and real part of this 
function for e = 1/5. Since Figs. 2 and 3 are drawn using a crude, uniform 
mesh of points for arguments and a uniform trapezoidal interpolation for 
shading of the function image, it happens from time to time that the mesh 
selects an argument point a bit away from a narrow peak (forming a wall 
of a priori infinite height) in the plotted function. The selected argument 
produces a small value of the function. When such small value is crudely 
connected with large values at neighboring points, a somewhat jagged shape 
is obtained. When the mesh spacing is decreased, drawings become smooth. 
However, the files with a small mesh spacing have a prohibitive size. Figs. 
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Fig. 3. The same non-local interaction density lunction as in Fig. 2, also multiplied 
by I/O -1 "! as in Fig. 2 to avoid suppression on the drawing when the distance between 
Xi and X2 increases, plotted for several choices of the points Xxi with X(\x^\,x~) — 
(1,6), (4,6), (6,4), and (6, 1), as indicated. See the text and Eqs. (43) to (48). 



1, 2, and 3 are provided as a result of a compromise between precision of 
rendering and size of the figure files. The compromise is made in such a 
way that in every case the key feature to be displayed and discussed is not 
altered in any significant way. 

3.3. Non-locality for A ~ m 

In attempts to understand mass generation in a relativistic theory of par- 
ticles, such as attempts to solve QCD for quark and gluon wave functions of 
hadrons in the Minkowski space or attempts to resolve the strong interac- 
tion structure of photons that results in mixing of photons with p-mesons, 
one is faced with a challenge of understanding dynamics of binding of effec- 
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Fig. 4. Modulus and real part of the function f(x) in Eq. (49). Shape of |/(x)| 
explains the somewhat jagged shapes of functions displayed in Figs. 2 and 3 ; see 
the text. 



tive constituents of some mass m. In QCD, the effective constituent mass 
m that dominates the mechanism of binding of lightest quarks appears to 
be of the order of Aqcd- The problem of chiral symmetry breaking in QCD 
and other theories can be rephrased as a question of what are the dominant 
interaction terms in an effective Hamiltonian of QCD when the decreasing 
RGPEP parameter A becomes comparable with Aqcd- In any case, unless 
the parameter A is comparable with or smaller than the constituent mass m, 
a large number of constituents may participate in the dynamics described 
by Hamiltonians evaluated using RGPEP. Therefore, it is of interest to con- 
struct and understand the structure of Hamiltonian interaction terms with 
small A, where by small it is meant that A is comparable with m. This sec- 
tion discusses the non-locality of first-order RGPEP three-particle vertex in 
this regime. 

It is visible in Fig. 1 that the RGPEP momentum-space form factor fx 
peaks at £ ~ 1/2, i.e., in the center of a middle bump in Fig. 1. Section 
3.1 shows that the two neighboring regions with ( below and above 1 
both contribute the same operator that the region < ( < 1 contributes. 
As A decreases, the size of fx also decreases. The smaller fx the weaker 
the interaction that changes the number of effective particles. However, 
in asymptotically free theories, the effective coupling constant gx increases 
when A decreases; gx may partly compensate for the small size of fx at 
its peak. The resulting interaction strength may thus be not eliminated 
entirely when A decreases to about m. Instead, the strength is located in 
the domain in momentum space where fx is maximal. Therefore, for A ~ m 
or smaller, given that the three regions of £ < 0, < ( < 1, and 1 < C 
contribute the same operator to the Hamiltonian, it is sufficient to consider 
the middle region of ( in Fig. 1, i.e., the region where fx forms a bump 
around z = ( ~ 1/2. In this region, the relative transverse momentum 
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q L = k 1 - is limited to values not exceeding order A, tempered in addition 
by how much the variable z deviates from 1/2. 

For a closer inspection of the middle bump region for A ~ m, consider 
the form factor 

with M.\ 2 given in Eq. (33), rewritten as 

f x = e -(3m 2 A 2 ) 2 e -[(A1? 2 -m 2 ) 2 -(3m 2 ) 2 ]/A^ ^ 

The exponential in front describes the size of the form factor at its maximum 
at z = 1/2 and q 1 - = 0, while the remaining factor describes the form factor 
fall-off away from its maximum. Now, 

(M 2 12 - m 2 ) 2 - (3m 2 ) 2 = (M 2 l2 + 2m 2 )(M 2 l2 - 4m 2 ) , (52) 

where the first factor on the right-hand side is not smaller than 6m 2 while 
the second factor can be small. These two factors can be analyzed in detail 
in terms of variables 

k L = q ± , (53) 
k z = (z-l/2)M 12 , (54) 

dz = 4z(l - z)4l~ , (55) 



M 



12 



with which 



z V ym 2 + k 2 / 
M\ 2 = 4(k 2 + m 2 ), (57) 



and 



,o 2/^2^2 £ 2 +3m 2 /2 k 2 
f x = e -(3m 2 /A 2 ) 2 g (A/2) 2 (W . (58) 

For A ~ m or smaller, the relative momentum of particles 1 and 2 is 
smaller than A/2 ~ m/2 and the RGPEP form factor can be very well 
approximated by Gaussian 

h = e - 9m 4 /A 4 e -24,n 2 fc 2 /A4 (59) 
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Using this approximation in Eq. (37) written in terms of variables R and r 
defined in Eqs. (B.l) to (B.4), one obtains 

^ I" 4z(l - Z)d 3 k e _24m 2 fc 2 /A 4 i[-{z-l/2)Pr+k ± r- L ] ^ 

J 2(2tt)3A^ 12 ' - y ' 

Since for A ~ m or smaller one has z ~ 1/2 and 

F = (z- l/2)Mi2 ~ (z- l/2)2m, (61) 
the function f\(xi, x 2 , x$) is approximated by 



f = 3e -9mVX* f d*PP+e{P+) iP{R _ X3) 

J 2(2vr)3 



f d 3 k 
J 4m(2vr 



_ e -24m 2 fc 2 /A 4 e i[k'(-Pr/(2m))+fc- L r- L ] _ ^ 



The key observation is that one can introduce a three-vector 

f = f(P, r ) = (^,r^ , (63) 



and write 



7 _ o -9m 4 /A 4 f d 3 PP + 6(P + ) iP r R _ X3 ) 



4m(2vr) 



2(2vr) s 

e -24m a W (64) 



so that the non-locality is parameterized in terms of r. Details of evaluation 
of the non-locality are described in Appendix C. Like in Eqs. (43) to (48) 
for massless particles, one can write 

fx(x 1 ,x 2 ,x 3 ) = X 6 g(p,T) , (65) 

where t 1 - = tp ± + sn 1 - and 

9(P,r) = S( 8 )^h p (t,T-). (66) 



Using 7 = Eq. (C.15) implies 

8 ft 1 /o\' 2 



i -5M*- 1 / 2 ) 
(r - tp -Me) IP I e 3^- 
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Fig. 5. Modulus of the non-local LF interaction density f\(xi,X2, £3) in Eq. (65) 
for massive particles with A = m, for which 7 = 1/6, using the same convention as 
in Fig. 2 but in terms of | p -1- 1 times h p in Eq. (67) multiplied in addition by e^ 1 /j. 
The extra multiplication is needed for removing the exponential suppression of the 
form factor h p so that the height of the plotted function matches with the one in 
Fig. 2. For further explanation, see Eqs. (65) to (67) and the text. 



For 7 = 1 (or A = \^6m), this result matches the massless case of Eq. (47) 
up to the factor e -1 / 4 when t ~ 1/2. In this case, t(l — t) amounts to 1/4. 

Non-locality of the interaction vertex for A = m is illustrated in Fig. 5, 
which should be compared with Fig. 2. For AX3 away from the points with 
t ~ 1/2, especially at the end-points, the result for A = m shown in Fig. 5 
is quite different from the massless case shown in Fig. 2; there is a strong 
suppression instead of enhancement at the end points. The smaller 7, i.e., 
the smaller A/m, the narrower the peak at t = 1/2 and lesser exponential 
suppression of the non-local vertex due to increase of \p^\. 
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In contrast to Fig. 2, Fig. 5 exhibits the feature that the non-local 
interaction density prefers the position of X3 in the middle of X\ and X2- This 
is a characteristic behavior for non-relativistic interactions, i.e., interactions 
in which relative momentum of interacting particles is smaller than their 
masses. Interactions of relativistic particles do not have this feature. 

Finally, one can observe that the three points xi, X2, and X3 must ap- 
proximately lie on a straight line on the LF. This feature is the same as in 
Fig. 3 and does not require a separate drawing. 



There is an analogy with a non-relativistic quantum mechanics that 
may be helpful in interpretation of the results described in the previous 
sections. The difficulty of interpretation may be expected because RGPEP 
is developed in the front form of Hamiltonian dynamics and the non-locality 
of interaction Hamiltonian densities is obtained on the LF hyper-plane in 
space-time rather than in space at a single moment of time. Thus, the 
relativistic RGPEP involves concepts that are general enough for hoping 
that the method will apply in derivation of wave functions of bound states 
of partons in theories with asymptotic freedom and infrared slavery or in 
attempts to understand symmetry breaking and mass generation at high 
energies. The space-time concepts do not appear quite intuitive from the 
point of view based on non-relativistic quantum mechanics. The latter only 
applies when motion of charged particles is characterized by velocities v ~ 
ocqedc, where c is the speed of light and olqed is the fine structure constant 
~ 1/ 137. This constant determines the strength of interactions which govern 
behavior of electrons bound in atoms. Binding of quarks and gluons occurs 
in space-time as a result of interactions with coupling constants about 100 
times greater than in QED and the non-relativistic intuition is not directly 
applicable. The same difficulty with building a physical picture is faced in 
all theories of mass generation and related symmetry breaking. 

Fortunately, an intuitive picture to think about can be arrived at by 
comparing RGPEP form factors f\ for small A with simple form factors 
that may be introduced ad hoc in an effective non-relativistic theory. For 
example, consider interactions that resemble emission and absorption of 
mesons by nucleons in nuclear physics. Suppose an effective interaction 
Hamiltonian has the form 



3.4- Non-locality in a non-relativistic theory 




(27TYS%p 1 +p2-p 3 ) 



X /A K Pl a \ P2 «AP3 + h - c -) > 



(68) 
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where pi,i = 1, 2, 3 are standard three-dimensional momentum variables 
conjugated with standard rectilinear co-ordinates x*i,i = 1, 2, 3 in space, 
respectively. Let an arbitrarily chosen form factor be 



fx = e -(A'P2) 2 /x\ (6Q) 

Suppose that one builds effective fields at the moment t = using opera- 
tors cl\ (only annihilation operators) and evaluates a non-local Hamiltonian 
density in space. Instead of Eq. (37), one is led to consider an expression 
of the form 

I, \ f d 3 P f d 3 q , _ lX 

fx( Xl ,x 2 ,x 3 ) ~ J j^J p/ A e , 

(70) 

X = (P/2 + q) Xl + (P/2 - q)x 2 - Px 3 (71) 
= P[( Xl + x 2 )/2 - x 3 ) + q( Xl - x 2 ) , (72) 

and obtains 

/a(xi,X 2 ,X 3 ) ~ S 3 [X 3 - (x 1 + X 2 )/2] A 3 e -A 2 (5i-5 2 ) 2 /16. (73) 

This result has an intuitive interpretation. Particle 3 is annihilated 
(created) exactly in the middle of positions x\ and x 2 where particles 1 and 
2 are created (annihilated), respectively. The width of the distribution of 
distances between points x\ and x 2 is 1/A. When A — > 00, the interaction is 
local. When A becomes small, it is greatly delocalized in the sense that the 
distance between x\ and x 2 can be large, but at the same time X3 is always 
precisely in the middle between x\ and x 2 . 

The situation is somewhat similar to the one shown in Fig. 5, except that 
in Fig. 5 one of the directions is on the LF. However, when A is small in 
comparison with masses and x 3 is fairly located on the LF between x\ and 
x 2 , one can observe that for all co-ordinates involved x + = and x~ = — 2x 3 . 
Therefore, in the non-relativistic case, r~ can be seen as analogous to — 2r z . 
This way one recovers the simple non-relativistic interpretation of the LF 
non-locality. Further discussion is provided in the next section. 



3.5. Relationship to a 2-body wave function 

Previous section provided an interpretation of the non-local interaction 
Hamiltonian densities on the LF by analogy that was limited to slow par- 
ticles. This section provides another intuitive picture that is more precise 
and not limited to slowly moving bound states. In fact, it can be used in 
any frame one wishes to use, including the infinite momentum frame (IMF) 
and the center of mass frame (CMF) alike. 
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In order to relate the non-locality of an effective Hamiltonian interaction 
term to a wave function of a 2-body bound state on the LF hyperplane 
defined by condition x + = 0, we introduce three species of particles which 
are annihilated by operators a±, a 2 , and 03, respectively. The three species 
are introduced to avoid the need for symmetrization of relevant functions 
for identical particles. Masses of the three species are assumed all equal m 
in order to match the conditions set in the previous sections. 

Consider the matrix element 

Mx u x 2 ) = (0\iP+ x ( Xl )iP+ x (x 2 )H XI al xp \0), (74) 

where |0) denotes the vacuum state that is annihilated by a±x p , a 2 \ p , and 
a^x-p for p + > 0. The effective quantum fields ipi\(x), ip2\(x), and ?p3\(x) 
are defined by Eqs. (19) to (22) for species 1, 2, and 3, respectively. The 
superscript + means differentiation, 

ij)+{x) = id + ip{x). (75) 

The state \P) with P + > 0, defined as a particle of the third species 

\P) = 4 Ap |0) , (76) 

can be considered analogous to a bound state of two effective constituents of 
species 1 and 2 corresponding to scale A. Namely, the eigenvalue equation 
for a bound state with momentum P, 

(H X0 + H XI )\P) = P~\P), (77) 

can be rewritten as 

\P) = p_ 1 n H XI \P). (78) 

The Hamiltonian H\j couples a pair of particles 1 and 2 with a particle 
3. Therefore, the eigenvalue equation for the state \P) involves both the 
2-particle component of type 1 and 2 and the 1-particle component of type 
3. This LF situation resembles the Lee model [16, 17]. The wave function to 
focus on is the 2-particle component. One can insert the identity in 2-body 
space, 

I12 = jd 3 x 1 d 3 x 2 ViA(^i)^ 2 A(a: 2 )|0)(0^+(xi)V 2 + A (x 2 ), (79) 
on the left-hand side of H\j and evaluate the matrix element 

Mx 1 ,x 2 ) = (OlV&foi) ^ A (s 2 ) p _ 1 rr I12 H XI \P) . (80) 
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The result takes the form 

tpxp(x 1 ,x 2 ) = J d 3 x' l d 3 x' 2 G\p(x 1 ,X2;x' l ,x 2 )<i)\p(x' 1 ,x 2 ) , (81) 

in which the matrix element 

Gxp(xi,x 2 ;x' l ,x' 2 ) = (0|^+ A (ai) i/>£ x (x 2 ) p ^ j VuK) V>2a(4) |0) , 

(82) 

appears in the role of a two-body Green's function for states with LF energy 
P~ on the LF hyperplane x + = 0. The matrix element of Eq. (74) plays 
the role of a vertex function in Eq. (81). 

The Hamiltonian H\j from Eq. (34), plainly modified to the case of three 
species of particles by associating position x& with species k for k = 1,2,3 
and keeping the factor 3 in front, can be inserted into Eq. (74), which yields 

M^,) = |V /Ae -«, +[< ,-„P- jM 

(83) 

By comparison with Eq. (37), the following relation is uncovered between 
the non-local interaction Hamiltonian densities obtained in first-order RG- 
PEP in previous sections and a 2-body vertex function on the LF: 

f°° dP + f d 2 P- L 
g\f\(xi,x 2 ,x 3 ) = J -^ 2 -<p\p(x 1 ,x 2 )e lPx:i . (84) 

Inverting the Fourier transform, one obtains 

<P\p(xi,x 2 ) = 9\ j d 3 x 3 f x (xi,x 2 ,x 3 )e- lPx3 . (85) 

Using expressions derived for f\(xi,x 2 ,x 3 ) in previous sections, one obtains 
in the case m = 

Mxi,x 2 ) = 3g x p+ e -iPR 

x [ 1 dzz(l-z)e-^- 1 ^ Pr -^ 1 -^ x2r±2 , (86) 
Jo 

and in the case A < m 

Mxi,x 2 ) = 3g x (J-^ P + e~ iPR 

x C(X/m) e "^[^^ +r±2 ], (87) 
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where 



C(A/m) 



e -9m 4 /A 4 



A 4 




(88) 



96 m 4 



Irrespective of the values of particle masses, RGPEP scale A, and the 
coupling constant g\, the vertex function contains the factor P + A 2 which 
carries its dimension. The vertex function contains as a factor a plane- wave 
function of the center-of-mass position variable R = {x\ + x-i)j2. Particles 
1 and 2 contribute equally to R because mi = = m. The remaining 
factor is a function of relative motion but depends on the total momentum 
P. More precisely, it has a universal property of being a function of two 
variables: square of r L = — and Pr. Taking into account that r + = 0, 
the former variable is equal to a square of a four- vector, — r 2 , and the latter 
variable is a product of two four-vectors. Both variables are invariant under 
7 kinematical LF symmetries. 

The generic structure of the vertex function implies that the Soper vari- 
able R s = z%\ + (1 — z)x2 [18] does not properly separate the center-of-mass 
motion from the relative motion of constituents. Namely, the parameteri- 
zation pi = zP + k and p2 = (1 — z)P — k implies p\X\ +P2X2 = PRs + kr, 
where the relative momentum k = (1 — z)p\ — zp2 has k + = 0. This means 
that the product PR S is a mixture of PR and Pr that depends on z. Pr and 
r 1 - can be arguments of the relative-motion vertex function, but z is not al- 
lowed to appear in the plane wave that describes the center-of-mass motion 
with a definite total momentum P. The product Pr is a natural variable 
to complement r 1 - as an argument of the vertex function even though these 
variables have different dimensions (see below). 

Another reason for the found structure of the vertex function to be of 
interest is that there exists an analogy between the LF wave function and 
AdS/CFT descriptions of bound-state form factors, discovered by Brodsky 
and de Teramond [19]. In that analogy, the effective transverse distance 
variable = \J z(l — z) r 1 - plays a key role as an argument of a bound- 
state wave function. The question concerning A-dependent non-local LF 
interaction Hamiltonian densities is whether their RGPEP evolution can 
be understood as dependence on a 5th dimension [20] in the context of 
duality [21, 22, 23] and whether this dependence can explain the shape 
of wave functions found by Brodsky and de Teramond. They interpret the 
analogy between AdS/CFT duality and AdS/QCD picture of hadrons in LF 
formulation without any need for considering the argument Pr and scale A. 

The issue of different dimensions of Pr and r 2 is resolved in the large-A 
case, i.e., A S> m — > and Eq. (86), by multiplication of r 1 - by A. In the 
small- A case, i.e., A < m and Eq. (87), the same issue is resolved by dividing 
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Pr by the sum of masses of the interacting constituents. The ratio of A to 
the masses becomes a dimensionless parameter. 

It is worth noting that the relativistic Eq. (87) predicts in the CMS a 
small difference between how = (r x , r y ) and r z enter the vertex function. 
Namely, for P 1 - = 0, the ratio Pr/(2m) becomes P + r~ /(Am) and P + in 
the CMS equals mass, say M, of the state under consideration (consider 
the Lee model [16]), while r~ = — 2r 3 when r + = 0. So, [Pr/(2m)] 2 in 
Eq. (87) becomes [Mr z /(2m)] 2 , or [1 - E B /(2m)] 2 r 2 , where E B denotes 
binding energy. This is how the LF vertex function keeps track of the mass 
defect due to binding. This result requires better understanding than the 
one offered here. On the other hand, the same Eq. (87) can be used in 
any frame, including the IMF, where one sees how the interaction vertex or 
wave function get squeezed due to motion. This is a relativistic squeezing 
in a quantum theory, not classical. 

Finally, it also seems worth mentioning that the relationship identified 
here between non-local LF interaction vertices and bound-state vertex func- 
tions may become helpful in calculating observables such as form factors. 
The suggestion is based on the fact that the function f\(x\, X2, £3) defined 
in Eq. (37) gives the same interaction Hamiltonian in Eq. (36) that is 
also obtained in Eq. (25) using f\(x\, X2, £3) defined in Eq. (26). Sup- 
pose that the old-fashioned perturbation theory for form factors [24] can be 
developed using non-local interaction Hamiltonians of the type defined in 
Eq. (25) with non-locality of the type defined in Eq. (26). Feynman rules 
with non-local vertices could then suggest how to incorporate all diagrams 
that count for all kinds of momentum transfers, not only those that have 
q + = 0. In the case of local vertex functions [25], it is known what to do 
when q + / 0. It is less clear what to do for bound-state vertices that in- 
volve non-trivial vertex functions. Therefore, it seems worth checking if the 
non-local interaction Hamiltonian densities obtained in first-order RGPEP 
lead to unique answers. 

4. Conclusion 

Renormalized Hamiltonian densities on the LF hyperplane in space-time 
contain interaction terms that are non-local. The non-locality is certainly 
intriguing and needs to be studied for many reasons as a feature of basic 
interactions. This article makes only a first step in this direction. In partic- 
ular, we calculate the non-locality in first-order of a perturbative expansion 
in powers of an effective coupling constant using RGPEP. This is done for 
terms that originate from a product of three fields. Such terms include cou- 
pling of fermions to gauge bosons, coupling of fermions to Yukawa particles, 
and coupling of non-Abelian gauge bosons with themselves. But the leading 
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non-locality is common to all these cases and can be calculated using scalar 
fields. The result is that the non-local interaction density has a generic 
form as a function of the space-time positions of effective particles that are 
created and annihilated by the interaction. This form can be understood in 
terms of a vertex function for a two-body bound state. 

The characteristic dependence of the vertex function on Pr and r 2 , where 
P denotes the bound-state total momentum and r denotes the relative po- 
sition of its two constituents, implies a characteristic dependence of the 
non-local Hamiltonian density on the position X3 of an annihilated (cre- 
ated) particle relatively to the positions x\ and X2 of created (annihilated) 
particles in a three-particle vertex. Namely, X3 is distributed near a straight 
line connecting x\ with X2 on the LF. Figs. 1 to 5 illustrate the shapes of 
the calculated distributions. 

All examples studied here were obtained for all particles having the 
same mass. For different masses, the results will be numerically different 
and of more interest from the point of view of application. However, there 
is no reason to expect a major alteration in the method and results beyond 
numerical changes. Nevertheless, such changes will be significant in practice. 

Another need for generalization concerns interactions that originate from 
products of more than three fields in one vertex. RGPEP provides a set of 
general rules for how to calculate non-local interaction densities in such cases 
order-by-order in perturbation theory. It is not excluded that the terms 
identified using RGPEP in perturbation theory will lead to a selection of 
dominant terms for which a non-perturbative evolution can be derived on 
a computer. 

A general speculation is in order concerning the role of non-local inter- 
actions in formation of strings of quantum gluons. Imagine that an effective 
gluon splits non-locally into two along a line as required by the LF non- 
locality of a Hamiltonian interaction term. Then, each of the two gluons 
interacts with neighboring gluons. It was argued before using analogy with 
RGPEP results for heavy quarkonia [26] that gluons of small A may attract 
each other in color singlets by potentials that resemble harmonic oscillator. 
In this context, the non-local splitting of gluons can be seen as a candidate 
for dynamical generation of quantum strings in which interactions between 
neighboring gluons in space are so strong that the string energy grows only 
linearly with its length, the effective gluon mass at scale A providing a unit 
of energy per unit of length order 1/A. 
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Appendix A 

Connection of Eqs. (36) and (37) with Eqs. (25) and (26) 

Details of the connection of interest are provided here for completeness. 
The connection involves several steps. Each of them involves manipulation 
of several variables. Understanding the connection requires tracing of these 
steps. These steps also exhibit an analogy between the momentum labeling 
of creation and annihilation operators in RGPEP and the parameters rf 
in the range between and 1 that appear in old-fashioned perturbation 
theory for scattering processes in the IMF [27], or parameters x in the same 
range from to 1 that appear in x + -ordered Feynman rules for calculating 
scattering amplitudes [14]. 

Creation and annihilation operators in the interaction Hamiltonian in 
Eq. (34) stand in normal order and inserting a sign of normal ordering does 
not change anything. Two terms, both with integration over P + > 0, can be 
changed to one term with integration over P + from — oo to +oo, rendering 
the result that H\j in Eq. (36) equals 




(A.l) 



One can use fx defined in Eq. (37) to introduce 



fx 



fx + fx ! 



(A.2) 



and write 




(A.3) 



where 




(A.4) 



X 



{zP + q)xi + [(1 - z)P - q]x 2 - Px z , 



(A.5) 




(A.6) 
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The function fx in Eq. (A. 3), denned by Eq. (A. 4), should be the same as 
the function fx in Eq. (25), defined by Eq. (26), i.e., 



fx{x 1 ,x 2 ,x 3 ) 



d 3 Pi 



JJy 2 (2tt)= 



2(2Tr) 6 S 6 ( Pl +p 2 +p 3 ) 



x g -(AA4 2 /A 2 ) 2 e +i(pix 1 +p 2 x 2 +P3x: i ) 



(A.7) 



In order to exhibit equivalence of fx in Eq. (A.7) and fx in Eq. (A. 4), we 
change integration variables in Eq. (A.7) according to Eqs. (31) and (32), 
from pi and p 2 to P12, ( and k, 



Pf = (Pl2 , Pi = (Pl2 + « , Pi 



Pi 2 + m 2 



CP, 



(Ai 



12 



pt = (i - qp+ , P i = (i - OAi - ^ , p 2 = ^ 2 _ + c) p+ • ( A - 9 ) 



rename p 3 to P, and obtain 



/ A (xi,x 2 ,x 3 ) 



/• rf 3 P r d 3 P 12 f +co dQP^l f 

J 2(2vr)3 J 2(2vr)3 2(2tt) J 



d(\P+\ f d 2 K A 



2(2*? J 2(2tt)3./_ 0O 2(2vr) J (2vr)2 
^3x3^D._ , d\ „-(AA1 2 /A 2 ) 2 -iY 



2(2tt) 3 5 a (Pi2 + P) e -l^"/*-r e -« , (A.10) 
^ = -(CP12 + «)ari - [(1 - C)Pi2 - k]x 2 - Px 3 , (A.ll) 

with AM 2 given by Eqs. (17) and (18). Change of the variable k 1 - to — k 1 - 
and integration over P\ 2 yields 



/ A (xi,x 2 ,x 3 ) 



d 3 P 
2(2^ 



'°° d(\P + \ f d 2 K^ (A>t3/A 2 ) 2 iy 
2(2vr) 7 (2tt)2 ' ' 



y = (CP + «)xi + [(i - C)p - k]x 2 - p X3 , 

which is to be compared with Eq. (A. 4) for fx(xi, x 2 , x 3 ) using 



(A.12) 
(A.13) 



p{ 



<p+ , pf = -CP X - ^ , P f 



p\ 2 + m 2 

-CP+ 



(A.14) 



pt = -(l - C)P + , Pt = -(i - C)P X + ^ , p 2 ~ 

p±2 + m 2 



pt 2 + m 2 



Ps 



P+ 



P3 



P3 



P+ 



-(i-C)P+' 

(A.15) 
(A.16) 
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and also Eqs. (17) and (18). 

Note that Eqs. (A. 12) and (A. 4) look similar in terms of integrations 
and functions they involve. For example, both involve the same integration 
over P and the functions X and Y coincide after replacement of z and q 1 - 
by C and k -1 , respectively. However, the two expressions as a whole differ 
by a factor of 3 in front. The integration over z ranges from to 1 while 
integration over £ ranges from — oo to +oo. This difference is analogous to 
the difference between parameters rj [27] or parameters x [14] that range 
from to 1 (rj and x were mentioned at the beginning of this Appendix) 
and the usual momentum variables that range from — oo to +oo. The form 
factor fx appears defined in terms of particle momenta differently in both 
cases. In order to see how these differences conspire to produce the same 
result, we first evaluate AM 2 in Eq. (A. 12). 

We divide the range of integration over £ in Eq. (A. 12) into three ranges; 
one from — oo to 0, shortly called 1, one from to 1, called 2, and one 
from 1 to oo, called 3. In each of these 3 regions, we change variables 
of integration from £ and to new variables and demonstrate that the 
resulting expression coincides with 1/3 of f\(xi,X2,X3) in Eq. (A. 4). 

Eqs. (17) and (18) imply 

AM 2 = (pi+P2+P3)(siPi + s 2 p 2 + s 3 p 3 ), (A.17) 
where Sj = sgn(p+). Using Eqs. (A. 14) to (A. 16), 

S3 = s P , s 2 = -SX-^Sp, si = S£Sp . (A.18) 

Thus, for all particles in the interaction term having the same mass m, 

AM 2 s P = m 2 {l - - si_ c ) 

+ (1 - s^pipz + (1 - si-c)P2P3 - {sq + Sl-<)Pl£>2 • (A.19) 



The resulting four-vector products read 



.20) 
.21) 
.22) 



PlP2 
PlP3 
P2P3 



k ±2 + m 2 



2C(1 " 



m 2 , 



( m z k ±2 + m 



2 2C ' 

(1 - C) m 2 k ±2 + m? 
2 2(1-0 ' 



(A 
(A 
(A 



7 , , f [ +0 ° d(\P + \ f d 2 n 



g -(A.M 2 A 2 ) 2 e -iY 



and the result is that 
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1M 2 ,2 K ±2 +m 2 ] 1 + ICI + I1-CI 



m 



C(i - 



2s P 



(A.23) 
(A.24) 



y = (CP + k)zi + [(1 - 0^ - k]x 2 - Px 3 , (A.25) 

should match the corresponding expression in Eq. (A. 4), repeated here for 
the readers' convenience, 



f\(x 1 ,x 2 ,x 3 ) 



d'P f 1 dz\P+\ /rfV (amVa2)2 , x 
(2vr)3 J Q 2(2vr) J (2ir) 2 



(A.26) 



- 1 - 2 + m 2 



AM' 2 = * - ' -m 2 , (A.27) 
z(l-z) 

X = (zP + q)x! + [(1 - z)P - q\x 2 - Px 3 . (A.28) 

The to 1 part of integration over £ in Eq. (A.23) matches exactly 1/3 of 
Eq. (A.26). The question is how to see the matching of the remaining 2/3 
with parts of integration over £ from — oo to and from 1 to +oo. 

There are two main regions of integration variables, P + > and P + < 0. 
Focus on the region P + > 0. Split integration over ( into three ranges; range 
1 from — oo to 0, range 2 from to 1, and range 3 from 1 to oo. In the 
region 2, results match. Consider region 1. For £ < 0, one has > 0, 
which means the particle 1 is annihilated, not created, particle 2 is created, 
and particle 3 is always annihilated for P + > 0. 

So, for C < 0, change variables treating p 2 as a total momentum (created, 
with a negative + component) composed of p\ and p 3 (both annihilated, 
with positive + components). This means 

p\ = -zp 2 -q, p 3 = -(1- z)p 2 + q, (A.29) 

and 

Pf = ~(P + = -zpt , Pi = ~CP ± - ^ = -zpi - q ± , (A.30) 

P + = -(i - ()p + -^-p + , pi = -(i- op x + ^ -> -p l , 

(A.31) 

p+ = P + = -(1 - z)p+ , pi = P ± = -(1 - z)p 2 L + q 1 - . (A.32) 
The required change of variables is 

C = j K_L = ~r~~ i (A-33) 



1 - z 1 - z 

dC = J T ^, ^ = ^ = (1-C)W. (A.34) 
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With this change of variables (sp = 1 here but it is kept as sp), 



AM Z 



m 



2 K ±2 + m 2 ] 1 + |C| + |1-C| 



C(i-C) 

q 1 - 2 + m 2 2 

~n r ~ m 



2s P 



1 



(A.35) 
(A.36) 



which is the expected function of z and q^. So, after the change of variables, 
the region 1 with P + > contributes 



2, n fd 3 P6{P+) f° d(\P+\ f d 2 ^ ~(am 



2 /a 2 ) 2 e -*y 



(A.37) 



d 3 P6(P+ ) f 1 dz\P+\ f d 2 q ± e -(AA4 2 /A 2 ) 2 _. 
2(2vr) 3 



?y 



2(2vr) 7 (2vr) 2 (1 - z) 4 



y 



l-z 

Change notation p 2 — > — P 

pi 



-zP + q P-q 

X\ + x 2 - Px 3 



- -(1-C)P + = -i* 



l-z 



pi 



(A.38) 
(A.39) 



-(1-QP ± + k ± 



keeping sp = sp, so that 



P = (l-z)P + q. 



(A.40) 



(A.41) 



Then, 



fi(x 1 ,x 2 ,x 3 ) = J 



d 3 P(l- z) 3 9(P+) f l dz\P+\{l 
2(27)3 

d 2 q ± e -(AA4 2 /A 2 ) 2 



2(2vr) 



J (2^)2 "713 



z) 4 



(A.42) 



where 



y = 



(-zP + g) xi + P x 2 - [(1 - z)P + q]x 3 . 



(A.43) 
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Changing variable P to — P, 

f( } = f d*P6(-P+) [Uj^P+l f dV e - ( AMW e -iY 

h(Xl > X2 > Xs > J 2(2vr)3 J 2(2vr) J (2vr)2 

(A.44) 
(A.45) 

Y = (zP + q) Xl + [(l- z)P- q]x 3 - Px 2 . (A.46) 

This result should be compared with 1/3 of f\{xi,x 2 ,x 3 ) given in Eq. 
(A.26). 

This comparison shows that initial integration range over P + > and 
thus particles 1 and 2 created and particle 3 annihilated in /, corresponds to 
integration over P + < and thus particles 1 and 3 annihilated and particle 
2 created in /. When both signs of P + are included in the integration, the 
result must be 

fi{x 1 ,x 2 ,x 3 ) = i/(xi,x 3 ,x 2 ) . (A.47) 

The next observation is that the integration over x 2 and X3 in the Hamilto- 
nian includes only the symmetric part of the functions f(xi,x 2 ,x 3 ) in Eqs. 
(25) and (A. 3). Therefore, Eq. (A.47) completes the explanation of how 
the integration over range 1 produces 1/3 of the interaction Hamiltonian. 

Since the integration in range 3 can be transformed in the same way, the 
only difference being that for £ > 1 the particle 1 and particle 2 are changed 
in their roles with respect to particle 3, it follows that the integration over 
( in range 3 produces the remaining 1/3 of the Hamiltonian. Hence, the 
connection of Eqs. (36) and (37) with Eqs. (25) and (26) is established, 
and both ways of writing the non-local interaction Hamiltonian that results 
from first-order RGPEP, are equivalent. 

Appendix B 

Integrals involved in non-locality for m = 
In terms of variables 

R = {x 1 +x 2 )/2, 
r = x\ — x 2 = x/X, 

which imply 

xi = R + r/2, 
x 2 = R-r/2, 



AM' 



q +m z o 
— m 



z(l-z) 



(B.l) 
(B.2) 



(B.3) 
(B.4) 
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Eq. (42) reads 



, d s P0(P + ) f 1 dzP + f d 2 q 



x e 

Integration over q 1 - renders 

72 „± D i2 



2(2vr) 3 7 2(2vr) 7 (2vr) 2 

P[.R+(z-l/2)r— rr 3 ]+i<2- L r- L 



7 (25?) 



A 2 z(l - z)tt (r -L A ^/I(T^)/2) 2 



r) 2 (2vr) 2 
and subsequent integration over P 1 - gives 

"°° dP+P+ f 1 A 2 z(l - z) 



(B.5) 



(B.6) 



/ A (xi,x 2 ,x 3 ) = 3 



47T 



47T 47T 



x 5 2 



R ± + (z-l/2)r ± -x~ 3 



x e 



~iP+[R- + (z-l/2)r--x-]/2-z(l-z)X 2 r ±2 /4 



(B.7) 



Integration over P + requires regularization. By inserting e eP+ / 2 under the 
integral and assuming that e — > 0+, one obtains 



/a(xi,x 2 ,x 3 ) 



-3A 



2 /•! 



cfe z(l — z) 



iGvr 3 7o [^7 + (1 



(B.8) 



x b 1 



ZXi + (1 — Z)%2 — X^ 



Hence, the function g(p,r) defined in Eq. (43) is given by Eq. (44). 

Appendix C 

Integrals involved in non-locality for A < m 

Introducing a dimensionless three- vectors p = Xf, where f is defined in 
Eq. (63), and p = k/X, where k is defined in Eqs. (53) and (54), one obtains 
the function / defined in Eq. (64) in the form 



f - . p -9myX* f d 3 PP+0(P+) lP(R _ 

1 ~ J 2(2vr)3 



X3) 



r x 3 d 3 P 

J 4m(2vr 



e a 



24m^ -2 
2- V & tPP 



(C.l) 
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Using /3 2 = 24m 2 /A 2 , one can perform integration over p, obtaining 

J 2(27r) 3 



A 3 v^ 3 



-p 2 /(2/3) 2 



/3 3 4m(2vr) 3 

To integrate over P -1- , / can be written in the form 
A 3 - 3 



(C.2) 



/ = 3e 



-9m 4 /A 4 



2 D_L 



d 2 P 

(2^ 



/3 3 4m(2vr) 3 



47T 



P+Ar~/2-P 



4m/3 



(C.3) 



Then, P ± can be written in terms of two mutually orthogonal transverse 
vectors, e p = p ± /\p ± \ and l ± , as 



P ± = (pet + qi^X. 



(C.4) 



The auxiliary parameter p has nothing to do with p = k/X introduced 
earlier. The integral over P 1 - becomes 



d 2 P A 



(2tt)* 



JP ± (R ± -xi) e 



X 2 5 



\i- L (R-x 3 y 



( P+\r~ /2-P ± p ± \ ' 
' I imp J 



2ir 



( P+Xr-/2-pX\p ± \ Y 
I 4m/3 J 



(C.5) 



.(C.6) 



The remaining integral over p, using x = e ~ p L (P ± ~ ^3")^) gives 



/ 



dp 
27 



e «pe^(R-a:3) ± A e 



P+Ar~/2-pA 



4m,9 



4m5 -(ZsexY yfa 
2| P -L| e Va| P -L|; A__ 

A|p-L| 2tt 



(C.7) 
(C.8) 



The entire integral is then 

A 4 



/ = 3e" 9m4 / A4 



(4vr/3)V 



2/3 



Ai- L (P-x 3 ) 



00 dP+ P+ e -i p+ (R--x 3 )/ 2 j£± 



47T 



(C.9) 
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The remaining integral over P + , after the same regularization by factor 
e _eP+ / 2 that was introduced in Appendix B, produces 



-9m 4 /A 4 \4 _f\p ± \ \ 2 _(m X \ 2 ( 2/3 \ : 

1 ~ 6 16tt3^| p ±| e 



K — Xn ±2 A r — te 



Xl- l (R-x 3 )- 



(C.10) 



Using dimensionless variables p = A (xi — x<i) and r = A (X3 — £2), one can 
write 



X = p ± (p ± /2-r ± )/\p ± \. 
Since t l must lie along p- 1 , so that t 1 - = tp ± , one gets 

X = -(t-l/2)|/|. 



(C.ll) 



(C.12) 



Thus, the approximate result for f\(x±, X2, X3) = A 6 g(p, r), where r J 
ip 1 - + sn^, is 

-1/4 



Wir 6 (t~ - tp- + ie) 



.£^_M!^ (t _ 1/2)2 9m4/A 4 



Writing 



one arrives at 



g( P , T ) = 5{s)r^h p {t,T ), 



1/4 e"^" — P-^- 1 / 2 ) 4e -9m 4 A 4 



(C.13) 



(C.14) 



/i p (t,r ) = 



(r- - tp- +ie) 2 
which results in Eq. (67). 



/3 2 



(C.15) 
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